Emergence of Symmetry in Complex Networks 
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Many real networks have been found to have a rich degree of symmetry, which is a very impor- 
tant structural property of complex network, yet has been rarely studied so far. And where does 
symmetry comes from has not been explained. To explore the mechanism underlying symmetry of 
the networks, we studied statistics of certain local symmetric motifs, such as symmetric bicliques 
and generalized symmetric bicliques, which contribute to local symmetry of networks. We found 
that symmetry of complex networks is a consequence of similar linkage pattern, which means that 
nodes with similar degree tend to share similar linkage targets. A improved version of BA model 
integrating similar linkage pattern successfully reproduces the symmetry of real networks, indicating 
that similar linkage pattern is the underlying ingredient that responsible for the emergence of the 
symmetry in complex networks. 

PACS numbers: 89.75.-k 89.75.Fb 05.40.-a 02.20.-a 



In the last decades, we have witnessed the great 
progress in the complex network researches [H, 0, [3, 0, 
i, m 0, S S, E, [ll E E, El- Previous studies 
have primarily focused on finding the statistical proper- 
ties of various networks, such as small world property [1, 
H) Hli power-law distribution of vertex degree build- 
ing block of network motifs^]; assortative mixing 



self-similarity [6[; community structure (10|, Il3l|and 
hierarchical structure of the network. And based 
on these properties, many network models, such as 
Barabasi- Albert (BA)[3 model. Watts- Strogatz modelQ 
have been proposed to help predict the future evolution 
of the network. However, an important property of net- 
work structure, symmetry, has been rarely studied. 

Concept of symmetry is based on the concept of au- 
tomorphism of the graph, which characterizes adjacency 
invariance to transformation operation on the node set. 
Graph has been widely used to represent systems con- 
sisting of components (represented by nodes) as well as 
their relation (represented by edges). If two nodes are 
connected by an edge they are defined as adjacent nodes. 
An automorphism acting on the node set can be viewed 
as a permutation of the nodes of the graph preserving 
the adjacency of the nodes. The set of automorphisms 
under the product of permutation forms a group |,2l] |. 
In general, a network is considered as asymmetric if its 
underlying graph contains only an identity permutation, 
otherwise, the network is symmetric. 

It has bee shown that various complex networks have a 
rich degree of symmetry 22, 2^. The fact that large real 



networks are symmetric is surprising [22l|. since 'almost 
all graphs are asymmetric' 26,] 28|. As an ubiquitous phe- 
nomenon, the existence of symmetry in the real networks 
strongly begs an explanation, since existing ingredients. 



such as continuous growth and preferential attachment^ 
dominating the construction of the network structure, are 
not dedicated to interpret the origination of symmetry in 
real networks. 

To explore the origin of symmetry in real networks, 
we summarize statistics of the local symmetric motifs 
contributing to the symmetry of the real networks, by 
which we found that similar linkage vattern f2^ 1. which 
means that nodes having similar property, for example 
degree, tend to have similar linkage targets, is a ubiqui- 
tous law that dominating the construction of structures 
of a variety of real networks. For example, in a friendship 
network, it is widely believed that persons with similar 
properties such as educational background, interest, age, 
would probably have common friends. 

To show that similar linkage pattern is a ubiquitous 
law that holds across many structures of real networks, 
we first summarize the statistics of symmetric bicliques 
in the real networks, which is a induced complete bipar- 
tite subgraph, denoted as Ky, , in which the degree of 
vertices in Vi is conserved [30|. Obviously, if graph G 
contains a symmetric biclique Kvi,V2i then the automor- 
phism group of G, denoted as Aut{G), will have a cor- 
responding geometric decomposition factor^^ Sn with 
n = \Vi\, which indicates that the size of Aut{G) has 
a factor of n\. Thus, Xvi^vi becomes a local symmetric 
motifsj23\ contributing to the symmetry of the network. 
Hence, symmetric bicliques will contribute to the symme- 
try of the network. Figure [T] illustrate two such bicliques. 
If we do not care about what Vi and V2 are, we also 

where \Vi\ = i and \Vj\ = j. 

Note 

that /Ci,i does not necessarily contribute to the local sym- 
metry of the network, hence, in the following discussion. 



use Kij to denote Kv^.V2j 



And the set consisting of all Kij is denoted as Idj 
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FIG. 1: Illustration of symmetric bicliques. Figure (a) shows 
an example of Kn,i, which contributes to the symmetry of the 
network with a geometric decomposition factor n\; figure (b) 
shows an example of K3^2, which contributes to the symmetry 
of the network with a geometric decomposition factor S3. 




Similar Patlern Size 
DR0(i=11 



Similar Pallern Size 
HOM(i=1) 



Similar Pattern Size 



FIG. 2: Size distribution of the symmetric bicliques for real 
networks. The horizontal axis for each panel is the size of 
symmetric bicliques and the vertical axis is the occurrence 
frequency of the symmetric bicliques with the corresponding 
size. Figure (a) and (b) show the biclique size distribution of 
the Internet at autonomous levelQ for /C„,i and /C„,2 respec- 
tively; Figure (c) shows the biclique size distribution of IC„.i 
of Homo Sapiens '2(§j. Figure (d) shows biclique size distribu- 
tion of /C„,i of Drosophila melanoaaster [20I] . 

"Here, the snapshot at 2006-07-10 of CAIDA[i1 is used. 



only JCn,i vi^ith n > 2 has been summarized. 

As shown in Table U , similar linkage pattern is a uni- 
verse phenomenon in the process of structure construc- 
tion of many real networks including social networks, bi- 
ological networks and technological networks. For in- 
stance, for ICn^i on InternetAS dataset, there are totally 
916 non-disjoint 31[ structures, among which there ex- 
ists some larger symmetric motifs, e.g. the maximal mo- 
tif has 343 nodes in Vi. For all the network we tested, 
simple symmetric motifs such as JCn,i and /C„,2 can be 
frequently observed. Moreover, for some networks, such 
as biogrid network 'DRO', even for some larger i, more 
complex symmetric motifs of JCn.i do exist. 

As shown in Figure [21 among those simple symmetric 
motifs with i = 1,2, the size (n) distributions are right- 
skewed with a long tail for larger size, which implies that 



a number of larger patterns do exist. 

Furthermore, we will show that similar linkage pattern 
will not happen in ER 27| random graphs. As shown 
in Table |TT1 only few randomized networks having the 
same size as the corresponding real networks, have sym- 
metric motifs of /C„,i; and no larger motifs /C„_i with 
i > 2 exist. Also we found that the number of motifs in 
JCn,i of randomized networks is much less than that of 
the corresponding real networks with the same size, the 
complexity of the motifs are much lower than that of the 
corresponding real networks. 

The frequent occurrence of complex Kn^i in real net- 
works and the unfrequcnt occurrence of complex K^^i in 
random networks strongly suggest that there exist some 
laws dominating the structure construction process of 
real networks. Consider the dynamic process of the net- 
work growth. We assume that at some time a new node 
V is added to the network, and a symmetric motif Ky^yz 
will arise. Thus, from the facts we have observed, it's 
reasonable to believe that v will attach to the existing 
nodes under the principle oi preferentially linking to those 
nodes to which other nodes in Vi attach. Since nodes in 
Vi have the same degree, it's reasonable to believe that 
nodes having the same degree will have the same linkage 
pattern. 

However, as shown in Figure [3l in real networks, nodes 
having the same degree tend to have only similar targets 
not exactly the same target sets, and these local motifs 
exhibiting non-exact similar linkage behavior also have 
chance to contribute to the symmetry of the network. 
Clearly, these local substructures are the generalization 
of symmetric bicliques, in the way that structures con- 
straint of the clique is relaxed from being complete bipar- 
tite to only satisfying that all the vertex of Vi have the 
same degree in the clique. Thus, this kind of generalized 
symmetric biclique can be denoted as Ky^ , where d is 
the degree of any vertex in Vi; in some cases that we do 
not care about V2, Ky^ or K^y^^ is often used. 

In a network, if non-exact similar linkage pattern does 
make sense, then nodes with the same degree will tend to 
share similar linkage targets. Thus, we need to measure 
in what extent these nodes share the same targets. Let 
V{m) = {v : V ^ V and d{v) = m} be all nodes with 
degree m, then the linkage targets of these nodes could 
be denoted as V{my = {v' : {v,v') G E and v e V{m)}. 
Then we could define 9m as EquationfTJ which is the ratio 
of the actual number of linkage targets of V{m) to the 
maximal probable number of linkage targets ( Maximal 
set of linkage targets could be obtained when overlapping 
of targets of nodes in V{m) is forbidden.) Obviously, 
this measure is a key index that can be used to quantify 
overlapping ratio of linkage target of nodes in V{m). 



\V'{m)\ 
m\V{m)\ 



(1) 



Obviously, we have < 0™ < 1. If is given. 
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TABLE I: Symmetric biclique statistics of a number of real networks. We measure the size of the networks by the number of 
nodes and edges, denoted by N and M, respectively. For each i < 7, the statistics of symmetric bicliques contained in /C„,iwith 
n > 2 is measured. We use a triple tuple {S, Min, Max) to show the statistics of lC„,i, where 5* is the number of non-equal 
substructures in ICn,i, and the Min, Max are the minimal and maximal size of symmetric bicliques, respectively. If /C„,i is 
empty, S — 0, and Min and Max is not available, denoted as For some larger i, we also enumerate corresponding statistics 

of ICn,i- 



/C„,iwith n > 2 



Network 


1 


2 


3 


4 5 


6 


7 


some larger i 


arXivflQ]'' 
Internet AS 


(135,2,7) 


(42,2,4) 


(17,2,3) 


(13,2,2) (11,2,2) (1,2 


2) (2,2,2) 


i = 16,(1,2,2) 


(916,2,343) (1057,2,285) (90,2,25) 


(9,2,4) (2,2,2) 


(0,- 


-) (0,-,-) 


(0,-,-) 


BioGridf20] 
















SAC 


(51,2,15) 


(7,2,5) 


(0,-r) 


(0,-r) (0,-,-) 


(0,- 


-) (0,-,-) 


(0,-,-) 


MUS 


(7,2,44) 


(8,2,12) 


(4,2,6) 


(2,2,2) (0,-,-) 


(1,2 


2) (0,-,-) 


(0,-,-) 


HOM 


(366,2,44) 


(53,2,12) 


(21,2,6) 


(5,2,2) (2,2,2) 


(1,2 


2) (0,-,-) 


i = 8, 10,21,(1,2,2) 


DRO 


(418,2,40) 


(16,2,11) 


(6,2,3) 


(6,2,3) (3,2,2) 


(0,- 


-) (3,2,3) 


i = 8, 10, 21,(2,2,2),i = 15, 27,(1,3,3) 
= 17, 18, 19, 23, 25 (1,2,2) i = 39,(1,6,6) 


CAE 


(245,2,47) 


(9,2,5) 


(1,2,2) 


(0,-,-) (0,-,-) 


(0,- 


-) (0,-,-) 


(0,-,-) 


USPowerGrid[8l 


(137,2,9) 


(25,2,3) 


(0,-,-) 


(1,2,2) (0,-,-) 


(0,- 


-) (0,-,-) 


(0,-,-) 



"Here, the snapshot at 2006-03 of HEP-TH (high energy physics 
theory) citation graph [Tgl is used. 
''Here, the snapshot at 2006-07-10 of CAIDA[i3| is used. 



TABLE II: Symmetry biclique statistics for a variety of ER 
random networks. For each network tested in Table HI we gen- 
erate the corresponding ER random networks with the same 
size using P4 Ji?Jf [23]. Two parameters are needed, the vertex 
number A'^ and the average degree z to ensure the generated 
random network having the same size with the correspond- 
ing real networks. Similar to the evaluation in Table U we 
also measure the statistics for each K.„,i, but only statistics of 
ICn^i is available for some random networks. Larger symmet- 
ric bicliques ICn,i with larger i, do not exist in these random 
networks. 



Network 


N 


z 


/C„ 1 with n > 2 


arXiv 


27770 


25.37 


(0,-,-) 


InternetAS 


22442 


4.06 


(62,2,3) 


BioGrid 








SAC 


5437 


26.86 


(0,-,-) 


MUS 


218 


3.65 


(1,2,2) 


HOM 


7522 


5.32 


(3,2,2) 


DRO 


7528 


6.69 


(0,-,-) 


CAE 


2780 


3.13 


(21,2,3) 


USPowerGrid 


4941 


1.49 


(231,2,3) 



we have 



\Vim)\ 



< 6m < 1. Note that, the lower 



IS, 



the more frequently similar linkage pattern will happen, 
while for 9m' & closer to 1, the exact opposite is true. As 
shown in Figure [H for small values of degrees, all tested 
networks tend to have relative small 9m, which strongly 
suggest that for these real networks, in the process of 
network growth, nodes with the same small degree tend 
to have the similar linkage behavior. 

It has been shown in the BAQ model that many real 
networks have power law degree distribution, which can 
be attributed to two basic ingredients: (1) Growth and 
{2) Preferential attachment. In BA model, new nodes will 




FIG. 3: Illustration of non-exact similar linkage pattern. All 
nodes in Vi, i.e., dark nodes, have the same number of link- 
ages. However, the linkage target of these nodes are not ex- 
actly the same. Motifs induced by V\ as well as their linkage 
targets, do not necessarily contribute to the local symme- 
try of networks. In Figure 3(a) the subgraph induced by 
V — {ni, He} will not lead to any automorphism, while 
the induced subgraph in Figure 3(b) will result in an au- 
tomorphism p = {vi,V4){v5,vq) and the subgraph induced 
hy V = {vi, ...,vs} in Figure 3(c) also contributes an auto- 



morphism p = {vi, V4){v5 , V6)(v7 , vs) to the symmetry of the 
graph. 



be continuously added to the existing networks, and at 
each time step, a new node is preferentially attached 
to fixed number of m (this number m is referred to as 
initial degree of the newly added node) existing highly 
connected nodes. However, ingredients about symmetry 
have not been considered in BA model and other net- 
work generation models. To reproduce symmetric net- 
works with power law degree distribution, we propose a 
new network model incorporating similar linkage pattern 
into BA's two ingredients. For this purpose, two modifi- 
cations on BA model's two principles need to be carried 
out: 

1. Newly added nodes are linked to the existing nodes 
not only under the principle of preferential attach- 
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FIG. 4: Figure(a),(b),(c),(d) show 6m distribution of arXiv, 
Internet, Biogrid SAC and Biogrid DRO, respectively. 



merit, but also similar linkage pattern. The latter 
principle implies that newly added node with ini- 
tial degree m tends to link to the targets to which 
existing nodes with degree m in the network are 
linked. 

2. Initial degree m of newly added nodes follows a cer- 
tain distribution instead of being a constant value. 
In BA and other existing models, initial degree is 
constant, while in the following study, we will show 
that in many real networks, m follows a certain 
distribution. 

The probability, denoted by 11, that a new node with 
initial degree m will be connected to node Vi not only 
rely on degree ki of node Vi but also depends on whether 
Vi belongs to V/(m). To incorporate the ingredient of 
similar linkage pattern into the basic BA model, we need 
to increase ll{vi) for those Vi belonging to V/(r7i). Hence, 
we define parameter a to control the relative significance 
of similar linkage pattern in the formation of network 
structure. Note that for a given to, V/(to) is not neces- 
sarily to be non-empty, hence the probability 11 would be 
defined in two cases: when V/(to) = 0, then 11 should be 
defined as: 



where ki is the degree of vertex Vi] when ^/(to) ^ 
should be defined as : 



(2) 



n 



n(i;.) 



if V, i Vl{m). 



, where a G (0, 1]. 

At some time step t, we may have V^'(to) = and 
n would be reduced to pure preferential attachment in 



FIG. 5: Distribution of initial degree of arXiv dataset 



terms of the value of degree. It will happen frequently 
in the initial stage of the network growth in our model. 
Because the abundance of degree is limited in the ini- 
tial stage, whatever the seed network is. For example, if 
the seed network contains isolated vertices, only degree 
could be found, if the seed network is an regular network, 
such as a complete network, we could also find only one 
degree in the network. 

Equation [3] has only one parameter a to control the 
relative significance of pure preferential attachment and 
similar linkage pattern. It's clear that the larger a is, 
the less impact of similar linkage pattern on the network 
will be. When a = 1, the model is reduced to the pure 
preferential attachment according to vertex degree. 

Obviously,the basic assumption of BA model is that 
any node except for those seed nodes has the same ini- 
tial degree. However, for some networks, especially social 
networks and technique networks, whose historical data 
about initial degree of real networks is available, we can 
easily find that the initial degree of real networks may 
be far away from a fixed value or a value independent 
of degree. For example. Figure [5] shows the distribution 
of initial degree of a citation network constructed from 
arXiv data set. From this figure, we can see that for 
larger initial degree, the frequency follows power law dis- 
tribution rather than a fixed value. 

Assume that we grows the network in the way follow- 
ing the principle of preferentially attachment with simi- 
lar linkage pattern. If the initial degree is constant, then 
each time a new node was added to the network, fixed 
number (m) edges would be introduced into the network. 
Thus, the local symmetric motifs would concentrate on 
those subgraphs with structure closer to Kl^. If m is 
very larger than 1, it's contradictory to the above ob- 
served fact that the larger n is, the less frequently ivTf 
tends to occur. 

Hence, it is necessary to extend the initial degree 
from a fixed value to some given distribution. From this 
perspective, BA model would be considered as a special 
case of our model in the way that F{m) is specified as a 
constant value. 

The algorithm of the model incorporating the ingredi- 
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ent of similar linkage pattern is the following: 

(1) Growth: Starting from a small number(no) of iso- 
lated nodes, at every time step, we add a new node with 
m edges that link the new node to m different nodes 
already present in the system, where m follows a distri- 
bution F{m) and m < rh, where rh is the upper bound 
of the initial degree m. 

(2) Preferential attachment with similar linkage pat- 
tern: The probability 11 that a new node will be con- 
nected to node Vi is defined by Equation [2] and [3l 

The above improved model based on similar linkage 
pattern needs just three input parameters (no, F{m),a). 
For the notational convenience, the model is doted as 
SLP{no, F{m), a), where 'SLP' is the abbreviation of 
'Similar Linkage Pattern'. 

To test the effect of similar linkage pattern on sym- 
metry of the networks, we first give some measures of 
symmetry of networks. The degree of the symmetry of a 
graph usually could be quantified by aa = \Aut(G)\^^. 
i.e., the size of the automorphism group. In order to com- 
pare the symmetry of networks with different sizes, /3g 
has been used to measure the symmetry in (2^ , which is 
defined as: 



(4) 



, where N is the number nodes in the network. /3g mea- 
sures the symmetry relative to maximal possible auto- 
morphism group of a graph with N nodes. Another sym- 
metry measure 7g is also given, which is the ratio of 
number of nodes in all those nontrivial orbits [s^fa set 
of equivalents nodes under automorphism operation) to 
the number of all nodes in the network. Specifically, let 
V = {Vi, V2, Vfc} be the automorphism partition un- 
der the action of {Aut{G),V) on node set V, 70 could 
be defined as : 



7G 



El 



<i<fcjyi|>i 



N 



(5) 



As shown in Figure 6(a)[ with a varying from 1 to 0.1, 
i.e., with more influence resulting from similar linkage 
pattern exserted on the network construction, the auto- 
morphism group size of networks increase several hun- 
dreds of orders of magnitudes. The inset (a) and (b) of 
Figure 6(a) also show that another two symmetry indices 
Pg and increase with the decrease of a. Such facts 
also can be observed from the Figure |6(b)|6(c)| and[6(d7 



Hence, it's reasonable to believe that similar linkage pat- 
tern is responsible for the emergence of the symmetry in 
the networks. 

If we remove the ingredient of similar linkage pattern, 
we will find that ingredient of preferential attachment 
with initial degree following a distribution will not nec- 
essarily reproduce symmetry of networks. 

As shown in Figure [7l when average degree (fc)[33] is 
small (close to 1), the network have higher degree of sym- 
metry. Note that those networks with (fc) closer to 1 tend 




(a) 



(b) 





(c) 



(d) 



FIG. 6: The effect of a and size on symmetry of networks gen- 
erated by the SLP model, (a) Simulation of SLP model with 
a varying from 0.1 to 1 in increment of 0.1. The horizontal 
axis of the figure as well as two insets (a) and (b) is a, the 
vertical axis of the figure is IgQc, the vertical axis of the inset 
figure (a) and (b) are /3g and 70, respectively. In the simu- 
lations, we use no — 10, m — 10, t — 10000 and employ two 
kinds of initial degree distributions. Blue square(n) shows 
an exponential distribution P{m) = a'y~"^ with 7 = 8; red 
circle(o) shows a power law distribution P{m) — am^' with 
7 = — 1. Figure (b),(c) and (d) show the growth of symmetry 
indices including IgacPc and 7g(%) of networks generated 
by SLP model. In the simulations shown in Figure (b),(c) 
and (d), a power law initial degree distribution is employed 
with m = 10 and 7 = —1, we fix no as 10 and vary a from 
0.1 to 1( The arrow shows the direction of increasing a). We 
vary t from to 5000 and capture the snapshots of the net- 
work every 50 units of time, thus we could get 100 samples 
of networks with linearly increasing sizes. Clearly, for all a, 
the growth of automorphism group size ac shows a obvious 
exponential trend; and decrease of (5g shows an power law 
trend; and limiv^oo Pg = o"(c^) .When a varies from 0.1 to 1, 
all three symmetry indices decrease. The inset of Figure (c) 
shows the amplified local plot, which clearly shows the fact 
that Pa decrease with the growth of a. 



to have the structure of tree and it is desirable that tree 
tends to have higher degree of symmetry. Such result 
conforms to the result reported in [2^ that BA random 
trees and uniform random trees have higher degree of 
symmetry. 

As shown in Figure 7(a)|7(b)|7(c) when (fc) increase, 
the symmetry of the network rapidly decays to a constant 
level t(7), which is determined by the slope of the power 
law distribution. Obviously, the steeper the initial degree 
distribution is, the higher the symmetry level is. When 
7 = 0, i.e., the slope of the double log distribution plot 
is zero, then symmetry of the network rapidly decays to 
zero or value close to as (fc) increases. However, with 
the slope becoming steeper, the symmetry of the network 
rapidly decays to an approximately constant value t(7) 
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(a) 



(b) 




TABLE III: Statistics of Kn,i in some SLP networks with 
power law initial degree distribution. All the parameters are 
set the same as Figure 7(d) In this table, we records the 
number, the minimal, the maximal size of the maximal size 
of ICn,i with 7 as one of {0,-0.5,-1,-1.5,-2} and (fc) as one of 
{3.5,4,4.5,5}. 



(k) 





-0.5 


7 
-1 


-1.5 


-2 


3.5 


(35,2,4) 


(76,2,11) 


(137,2,10) 


(222,2,14) 


(295,2,26) 


4 


(35,2,5) 


(85,2,5) 


(149,2,6) 


(240,2,9) 


(281,2,41) 


4.5 


(36,2,4) 


(67,2,10) 


(121,2,9) 


(192,2,18) 


(251,2,47) 


5 


(26,2,4) 


(57,2,5) 


(113,2,6) 


(220,2,16) 


(246,2,40) 



(c) (d) 

FIG. 7: The effect of average degree (fc) on symmetry of net- 
works generated by preferential attachment with initial degree 
following a power law distribution P{m) — am'' . In this sim- 
ulation, we use a = 1 to eliminate the influence of similar 
linkage pattern, and other parameters are set as no — 10, 
t = 5000, and 7 = 0, —0.5, —1, —1.5, —2. For each 7, we vary 
(fc) from 1 to 5 in increments of 0.5. Figure (a),(b),(c) show 
the trend of automorphism group size Igaa, Pa and 7g(%) 
with the growth of the average degree of the network, respec- 
tively. It is clear that symmetry of the network will rapidly 
(super linearly) decrease to an constant level for less steeper 
initial degree distribution. Figure (d) shows the relation be- 
tween the number of local structure ifn,! and the slop (I7I) of 
the power law initial degree distribution for average degree as 
one of {5,4.5,4,3.5}. Parameters in Figure (d) are the same as 
Figure (a),(b) and (c). Clearly, with the increase of |7j, more 
Kn^i will occur as the substructures of the network. 

that is far larger than as (fc) increases. Thus, for steep 
log-log initial degree distribution, non-ignored degree of 
symmetry would be observed. Note that steeper initial 
degree distribution will result in a higher probability of 
smaller initial degree m, especially to = 1; as a result, 
more tree-like symmetry will be found in the structure 
of the network. As shown in Figure [7 (d)[ the number of 
Kns increases with the growth of I7I; as observed in Ta- 
ble [TITl the complexity and the size of Kn,i also increase 
with the slope of the double log initial degree distribu- 
tion. 

Thus, it's rational to conclude that only preferential at- 
tachment with initial degree following a distribution, will 
not necessarily reproduce symmetry of networks. Only in 
those cases that small initial degrees have higher proba- 
bility of occurrences, especially m — 1, will produce tree- 
like symmetry of networks. To reproduce higher proba- 
bility of smaller initial degree, we need to decrease the 
maximal initial degree or increase the slope of the initial 
degree distribution. 

In summary, through the statistics of certain local sym- 
metric motifs, such as (generalized) symmetric bicliques 



in many real networks, we found that similar linkage pat- 
tern plays an important role in the origin of symmetry of 
networks. To incorporate this ingredient into BA model, 
we improved BA model in two aspects: (1) extending 
the initial degree from a constant value to a distribu- 
tion; (2) increasing the linkage probability of those target 
nodes. Simulation shows that similar linkage pattern is 
responsible for the emergence of symmetry of networks, 
while preferential attachment with initial degree follow- 
ing a distribution will only reproduce tree-like symmetry 
in some cases. 

Extensive existence of similar linkage pattern in real 
networks inspires us that behavior of individual nodes is 
far away from randomness, which provide us a brand 
new perspective, symmetry, to imderstand the self- 
organization of the complex systems. From this new 
viewpoint, in our studies we have found strong positive 
correlation between similar linkage pattern and symme- 
try of networks, which demonstrates that the emergence 
of complexity at macro system level is originated from 
the simple micro mechanism of individual component of 
the network system. Since symmetry breaking is the ba- 
sic mechanism underlying the procedure of the network 
growth, we believe our studies will be of great value to 
help explore the laws dominating symmetry breaking in 
complex networks. 
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